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Abstract 

A diffuse matter filled Type V Universe is studied. The anisotropic behaviour, 
the distortion caused to the CMBR and the parameter region allowed by 

■ present cosmological bounds are examined. It is shown how the overall sky 
pattern of temperature anisotropics changes under a non-infinitesimal spatial 

■ coordinate transformation that preserves the Type V manifest homogeneity. 
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I. INTRODUCTION 



The general solution of Type V spatially homogenous cosmology in the realistic case of 
untilted diffuse mater p = g p, [|l]], 0, is studied in various forms of the metric. This 
kind of equation of state covers the actual matter composition of the Universe for most of 
it's time evolution. It covers both the radiation era and the dust epoch. Note also, that 
not only theoretical considerations (Mach's principle) but also observational data, suggest a 
value close to zero for the rotation of the Universe M , thus validating our assumption about 
untilted fluid. This solution describes an open Universe, a feature that is lately favorable 
H, and at the same time is not trivially simple like the Kasner family. 

The existence of a primordial global anisotropy in our Universe is an open issue. The 
promising string theories suggest an anisotropic expansion that compactified all but four 
dimensions. It is natural to expect that the remaining dimensions unfold to one of the nine 
known types of spatially homogenous models, most of which are anisotropic. 

On the other hand it is well known that the smoothness of the cosmic microwave back- 
ground radiation (CMBR) suggests that the Friedmann- Robert son- Walker (FRW) cosmo- 
logical model is a fair description of the actual present Universe. However the radiation 
spectra has small fluctuations originating, in general, from pertubations produced by some 
physical mechanism (inflation-topological defects) and from primordial anisotropies and in- 
homogeneities. If the Universe experienced an inflationary period, any primordial overall 
anisotropy may have been smoothed but not erased. A possible realistic scenario then, 
could be that we had, after the Planck regime, an anisotropic Universe that appeared more 
isotropic after inflation and almost FRW during the last scattering and afterwards. It is 
certainly important to know more about the characteristics and the evolution of this global 
geometric primordial anisotropy. 

The above mentioned scenario can be fully studied, only if we study the general solutions 
of the nine types of spatially homogeneous cosmologies. Interesting works have appeared 
in the literature regarding the connection of the various Bianchi types with the observed 
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cosmological data. There have been works assuming small anisotropic pertubations to the 
FRW solution [||], 0, [[7] which are valid after last scattering and others ||, 0, which 
incorporate the exact solutions for various types of homogeneous models. All these works 
study the final temperature patterns induced on the CMBR. The basic conclusion is that, if 
someone wants the distortion of the null geodesies after last scattering to be small in order 
to comply with the observed quadrupole, then one has to fix the parameters of the various 
models to produce a metric very close to the FRW. Consequently, they set upper limits on 
the relative shear (§) , where £ is the shear anisotropy and 6 = 3H is the expansion rate, 
and on the vorticity. 

There is also the important work of W. R. Stoeger, R. Maartens, and G. F. R. Ellis 
1(J, |jTTJ in which a theorem is proven stating that if all observers measured CMBR to be 



almost isotropic, in our past light cone from the last scattering till now, then this universe 
has to be almost FRW in this region. They are referring to an expanding Universe with 
noninteracting matter and radiation. Then based on this theorem and some not easily 
testable simplifications they set upper limits on the allowed relative shear (§J Q • 

One can in principle use the CMBR data in order to determine the allowed parameter 
space and the type of the various homogenous cosmological solutions. Assuming that the 
actual Universe after last scattering is a dust Bianchi V one, we studied the distortion of 
the initial isotropic temperature pattern caused by the null geodesies. We present a method 
that yields simple system of equations for the null geodesies (a similar idea has been used in 
fl2|l). We found that for some parameters that make the solution close to the FRW one, it is 
possible to have for the period from last scattering till now very small quadrupole anisotropy 
compatible with the present data. For parameters not close to the FRW solution we always 
failed to find acceptably small quadrupole anisotropy. The most important feature we found 
is that a non-infinitesimal and not orthogonal transformation of coordinates that preserves 
the Bianchi V manifest homogeneity of the spacetime, leads to a non trivial transformation 
of the repeatedly appearing in the literature overall sky-pattern of anisotropies, although 
does not affect the covariant properties carried in that pattern. 



The paper is organized as follows: In section II the Einstein equations are setup and 
completely solved. We present the solution in the off diagonal form, in order to make 
manifest the dependency on the spatial coordinate system, of the calculations leading to 
the temperature fluctuations as a function of the sky angles. In section III the anisotropic 
behaviour and its effect on the temperature pattern, meaning the plot of the temperature 
as a function of the sky angles, is studied. The magnitude of the pattern distortion resulted 
by the use of the non-diagonal form of the metric is presented. Some concluding remarks 
are also given. 

II. THE UNTILTED DIFFUSE MATTER BIANCHI V COSMOLOGY 

Einstein equations are written as 



where k = —8tvG the Einstein constant and Tab the energy-momentum tensor. Capital 
Latin indices are spacetime indices A,B,.. = 0,...,3 while lower Latin are spatial indices 



In the 3+1 analysis, we use the lapse and shift functions in order to write the 4-metric 

as 



{a) Ra 



- (4) i? g A B = kT A b, 



(1) 



1,2,3. 



ds 2 = (NiN* - N 2 )dt 2 + 2N i dx i dt + g^dxK 



(2) 



Then, Einstein equations (H) are equivalent to the following system : 




K\i = -kNTI 



tO 



(3) 




(4) 



k) -NKK) + N ^R) + g il N bl - (K^N 1 + K\N\ S - K]N\^ 




(5) 



The extrinsic curvature is given by 
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1 

2iV 



(6) 



while 



K) = g a K lh K = K\ 



where g u is the inverse of the 3-metric g^. We also define 



(7) 



(3) R i = g il (3) Ri 



T A — n ACrp 

l B = 9 I CB, 



(4) T = 



(8) 



The class of spatially homogeneous spacetimes is characterised by the existence of an Tri- 
dimensional isometry group of motions, acting on each surface of simultaneity S t . If m > 3 
and there is no proper invariant subgroup of dimension 3, the spacetime is of the Kantowski- 
Sachs type and will not concern us further. When m = 3, the dimension of S t , there 
exist three basis one-forms a a satisfying 



do" = Aa^ - a% = 2Cfrr?of, 



(9) 



where (7g are the structure constants of the corresponding isometry group. Greek indices 
number the different basis 1-forms and take values in the range 1, ...,3 . In this case there 
are coordinates t, x % such that the line element 13 takes the form 



ds 2 



N a (t)N a (t) - N 2 (t) dt 2 + 2N a (t)o?(x)dx i dt + laP (t)a°(x)a^(x)dx i dx j . (10) 



Inserting the line element Eq. fllOD into the system of equations Eq.(]3|), Eq.(f|), Eq.(j5]) we 
get, via repeated use of Eq.@ the following set of ordinary differential equations for the 
Bianchi-type spatially homogeneous spacetimes: 

kN. 



rpO 

2 a ' 



(11) 



K^Kl - K 2 + (3) R + 2 A = -2k (t o ° - T° p N p 



(12) 



N 



Kp -NKK* + NR« + 2N" {k«C$ p - K$C« p ) = (t« + T°N a ) - - [k(T ° + If) 



6 P 



(13) 



where 



K% = j ap K pP , m = T°»R pl 3 and (14) 



K. 



a/3 



^ {lap +2lauC% p NP + 2 lfSv C» p N») , (15) 



2N 

1C x pR C%laxl KV + 2C x aK C» vlpxl ™. (16) 

The components of the energy- momentum tensor appearing in the system of Eq . (|TT|) , 
Eq. (|i~2"D , Eq. (|i"3|) , are the time-dependent parts of the corresponding full components Tg. 
The energy-momentum tensor of a perfect fluid can be written as 

T AB = (p + p)u A u B + pg AB , «V = -1, (17) 

where u A is its four- velocity {ua = 9abu b ), p is the total energy density and p the pressure. 
The conservation equations T AB -b = are equivalent to 

p+Q(p + p)=0 (18) 
(p+p) u A +p iB g AB +Pu A = 0, (19) 

where O = —K = u A A is the expansion scalar and the dot denotes \j u . Now the last 
equations and the system Eq. (|3|), Eq.([|), Eq.@ are integrable given an equation of state 
p = p(p). A very general equation of state is that of the diffuse matter 

v = gp, \g\ < i (20) 

with g a constant. It includes the dust case for g = 0, gravitational waves (g = 1/3), 
cosmic string networks (g = —1/3), domain walls (g = —2/3) and free massless scalar fields 

(<? = +!)■ 

In spatially homogeneous spacetimes, if the perfect fluid moves orthogonal to the spatial 
slices, we have 

(2i » 



in the so called sychronous coordinate frame (N a = 0) which in our case it is also commoving. 
Using now Eq.flTo]) we can write Eq.fllSD as follows : 

(22) 



P 1 

+ ^ = 0, 



p + p 27 

where 7 = det (jap) ■ The dot denotes ordinary differentiation with respect to the time label 
that characterises the slices S 4 . Integrating the above expression with the help of Eq. (|2~0|) 
we find 



g+i 

P = Pil 2 



(23) 



The quantity pi is a positive constant. Note that Eq.(|r9|) becomes an identity since p = p(t) 
and u A is a geodesic field. The fact that our perfect fluid is taken to be untilted leads to the 
simplification T° = and thus the linear equations Eq.flTTD coincide with the vacuum ones, 
i.e. they are (when N a vanishes) 



7^ l vP -y* C£, 7 a/3 = 0. 



We now specialise to the case of the Bianchi type V, i.e. we set C\ 



(24) 



13 



-a 



31 



a 
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— C| 2 = 1/2 , all other vanish. It can be shown with the aid of a time independent automor- 
phism inducing diffeomorphism (A.I.D.) [|14j], that we can choose, keeping the full generality, 
7i3 = 723 = 0, 7n722 ~ (712) 2 = (733) 2 - It can further be shown that, after completely in- 
tegrating Einstein equations at this stage, the off diagonal element can be rotated away 
through a transformation 7 y a/3 = A^A^ 7^ with A a constant automorphism matrix of the 



form 



r s 



A 



t u 
1 



(25) 



Such a transformation emerges as the effect on 7 a/ 3 of a time independent "frozen" A.I.D. 



see Appendix of [ 375 1 ; Thus the appearance of a non vanishing 7 12 is tantamount to the 
use of a spatial coordinate system different from the system supporting the diagonal form 



of the solutions. In what follows we keep the non diagonal element non zero in order to 
solve the geodesic equations in this different spatial coordinate choice and thus be able to 
discuss the dependency of the temperature pattern on this choice. One would wonder why 
we should make the endeavour to solve the Einstein equations in this off diagonal form and 
not just take the diagonal solution and obtain a non diagonal form through an arbitrary 
time dependent automorphism. The reason is that only the constant automorphism of the 
form discussed above, keep us within the space of solutions, as it can be proven using the 
off diagonal solution presented here. This off diagonal solution represents the full space of 
solutions (with 713 = 723 = of course). The metric 7^ is therefore taken to be 



a 



lap 



b c 
/ 



l = f^ac-b 2 = f 



(26) 



and Eqs. (|24|) are then identically satisfied. An off-diagonal matrix of the above form gives 
the following Ricci tensor Eq . (jlB]) and Ricci scalar : 



Rap = 27~ 1/3 7«/3, 



Mr = 67- 1 / 3 . 



(27) 



As far as time is concerned , we adopt the gauge fixing condition N = y/j , since this 
simplifies the form of the equations. Consequently, the remaining Einstein equations Eq. ([12]) , 
Eq.QI3D read 



■ ap- 1 

7 l a p + l- 



- 24 7 2/3 + 87KP = 

lap -Y" i al xipv +G(7)7a/3 = 



(28) 
(29) 



with 



G(j) = -4 7 2/3 + 7 k (p - p) 



(30) 



Taking the trace of Eq. ( p9|) , one arrives at 
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-I + 3G( 7 ) = 



(31) 



which has a first integral 



where 



w = (— ] + UOy) = const, 
\7/ 



U{i) = 6 / ^-dj. 
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We make the conformal transformation 



(32) 



(33) 



™a/3 = 7 1/3 7a/3 



and then, due to Eq.(p7[) and Eq.(|3T]) we can write Eq.(^9|) as 



ZV aj 3 —w^ v U7 afl U7p u — 0. 



The above equation can be integrated immediately , giving 



(34) 



(35) 



(36) 



Since w = det(zu a/ 3) — 1 =>■ = . The system of Eq.(^) is equivalent to the following 
set of equations 



(7 



-1/3, 



(7- 1/3 c)- 



^( 7 - 1 / 3 a)+^( 7 -V3 6) 

^( 7 -l/3 6)+1? 2 (7 -l/3 c) 

^( 7 - 1 / 3 a)+^( 7 - 1 /3 6) 



(37) 
(38) 
(39) 
(40) 



and also the relations ?9f = t? 3 , = $3 = = $1 = 0. Then $| = — $\ . At this stage we could 
invoke a constant automorphism matrix of the form earlier given and set the matrix dp to 
the form 



' <&\ 0^ 



1? 









(41) 
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which would directly led us to the diagonal form of the solution (T|], ||, [|] (or Joseph's |16| , 
in the vacuum case). However since our aim is to uncover the dependency of the temperature 
pattern on the spatial coordinate system, we choose to keep the non essential parameters in 
the game. The quadratic Einstein equation Eq.(|28|) and the system Eq. fl3~T|) - Eq.pO|) imply 

w 1 



{®\f + m = |- ±l/( 7 ) - 12 7 2 / 3 + 4 7K p, 
where the function Ufa) for the diffuse matter equation of state Eq.(^) is written as 



(42) 



Ufa) = -367 2/3 + 12/tp l7 ^ 



1-9 



(43) 



Then Eq. (|42|) reduces the independent constants : 

For b , 7^ 0, we have from Eq . (|38|) , Eq.(^) and Eq . ([26|) that 

§\a 2 -2d\ab-§\b 2 = d\ 1 2 l\ 



(44) 



(45) 



Since Eq.(f40j) can be derived from the other three of the system Eq . (p7|) - Eq.(|4~0[), it suffices 
to solve the system of Eq. (|37D, Eq. (|39|) and Eq. (|45|) . The solution of the system Eq. fl37[), 
Eq.(p9l) gives the expressions 7 _1 / 3 a , 7" 

-i/3 6 

as sums of two exponential functions of time. 
Only for w > these solutions satisfy also Eq.(^). Thus finally, the general solution of the 
Bianchi V universe with an orthogonal perfect fluid and linear equation of state, in the off 
diagonal form, is given as follows : 



ait) = 7 1/3 («) 

b(t) = 7 1/s («) 

c(t) = 7 '/ 3 (t) 

f(t) = 7 1/3 («) 



" + S eXP (-V3 * 



|(vf-0-(vf')-SI(vf-0-(-vl' 



fc-'( g ) 



/i(x) 



+ 36a; 4 - 12/cpi 



.7; 



(46) 



The solution contains five arbitrary constants k, $\, fi 2 ,^ , pi . The integration constant 
that arises in the expression giving 7 (t) is absorbable (non-essential) and is fixed by choosing 
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the zero point on the time axis. One can further realise that three of the five constants 
are also absorbable as explained above, leaving us with two essential parameters one for 
the matter and one for the gravity, as expected. The counting of the number of essential 
parameters for the vacuum Bianchi cosmologies can be understood as follows: From the 
initial twelve constants 7a/?(t),7 a /3 (t ) we subtract twice the number of independent first 
class constraints and the number of outer automorphism parameters. In our Type V case 
this reach 12 — 2x4 — 3 = 1 essential parameter. For Type II we have 12 — 2x3 — 4 = 2 
and likewise for other Types. 

If one employs the above derived solution to describe the dust filled Universe, after last 
scattering, it is possible to constraint some of the free parameters. Thus 

w = ^T3> ( 47 ) 



-1 Sv^fi 

Pl = i^7*' (48) 



l 

where q is the decceleration parameter, H is the Hubble rate and Q is the flatness parameter. 
In their definition the mean scale factor 7 , has been used. The above expressions will be 
used to calculate the constants w, pi, jtoday = 7o from the measurable H Q , q Q , Q Q . 

The above relations impose the following constraints on the allowed parameter space of 
the solution: 

3fi 

2-q- — >0, <?"2 >0 - ( 5 °) 
The above constraints imply that < q < 2 and < Q < 1. 



III. THE ANISOTROPIC BEHAVIOR 



It is well known flTTfl that, under the assumptions of zero cosmological constant, the 
validity of the dominant energy condition and the positiveness of pressure, Bianchi types I, 
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V and VII are the only ones that can satisfy four isotropisation conditions enlisted below. 
Of course this by no means implies that the solution will indeed isotropise. We have to know 
the exact solutions in order to examine the anisotropic behaviour. These four conditions of 
C. B. Collins, S. W. Hawking, are as follows : 

(i) 7 — > oo 

(ii) T 00 > and f£ -> 

(iii) § — > , where 2£ 2 = a AB a AB 

(iv) I' 1 ^ la/3 const 

The limits are understood to be taken as the proper time approach infinity r — ► +oo. 
Referring to the first condition let us study the behaviour of the quantity j(t), for the metric 
presented in the previous section. The proper time r is given, in the dust case and for an 
expanding universe, by the following relation 

r(,) = j[F(x)d», F (x )^ 7 ==^= i , (51) 

where x = 7 1//6 and the choice of the origin of time is r = for 7 = 0. It is easily verified that 
7 — > 00 for r — > +00. The second condition is trivially satisfied due to the orthogonality of 
the four velocity u A . For the third criterion, we can compute the shear anisotropy S via the 
shear tensor a AB in the basis where N a = 



2S 2 = K«Ki -lK* = -L: (A - MR + 2.T °) . (52) 



Finally, we find for the solution Eq.(|46| 



E 

e 



w 



1/2 

(53) 



_3w-U(-f) 

One should note that for 7 — >• we have ^ — > Furthermore, the function ^ is a strictly 
decreasing function of 7 and lim-^oo ^ = 0. The last condition (iv), holds for dust, since 
we can prove that lim x ^ +00 t(x) < +00 . Therefore the solution satisfies all the specified 
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conditions. In the following, we will examine some more criteria, that can characterise the 
behaviour of a homogeneous model. 

A cosmological model with a perfect fluid source is said, according to GFR Ellis |18 



to be close to a FRW model in some open set of the spacetime, if for some suitably small 
constant e C 1 the following inequalities hold, in this open set: 



where E\ v = Cx^vk rfrf, H\ v = Cx^uk ff'rf are the electric and magnetic parts of the Weyl 
tensor Cx^ UK respectively. The dual Weyl tensor Cx^uk is defined as Cx^uk = ^Vab^stcd with 
r) abcd the complete antisymmetric tensor and r] 12U = . These covariant criteria are 
believed to indicate that the model under consideration is close to a FRW one. For the 
solution considered here we find that 



E^E, V \ J\2w + 3U\ ^\H^H, U \ y/3 

6vj — ~ — 7t ' m = 6 vGw- — -. (55) 



H 2 w-U ' H 2 w-U 

These inequalities always hold in some interval (71,72) if we choose a sufficient small value 
for the parameter w. 

Finally the distortion that the geometric anisotropies induce to the temperature pattern 
via the null geodesies will be examined. Through this analysis it will be possible to find 
also the parameter region that is allowed from the current CMBR constraints. The small 
anisotropies of the CMBR are studied computing the angular variation of the microwave 
temperature, T o (0, (p) characterising the frequency distribution of the radiation. One defines 

where (T ) is the mean, over the sky, observed temperature. The above field is decomposed 
in spherical harmonics 

-y-(M = E E «i»*WM- (57) 

1=1 m=—l 

The various multipole moments ai m are defined by 
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r2n rn AT 

ai m = dip dd sin 9 -—Yi m (9,(p) (58) 
Jo Jo 1 

and we can finally define 

1 1 

(cn) 2 = (|a /m | 2 J> = — \ a im\ 2 ■ (59) 

m=—l 

After subtracting the dipole moment of anisotropy which is ascribed to the local motion of 
our sun and our galaxy around the local supercluster, one finds a small anisotropy of the 
order of 10~ 5 . Part of the remaining anisotropy is the quadrupole fluctuation. 

The cosmological redshift z e of a photon emitted at last scattering can be expressed as 



(u A k A ) 

1+ * = k*v (60) 

where k A is the tangent vector to the null geodesies and u A is the four velocity which in our 
case is vertical to the homogenous slices. The subscripts denote the emission time and the 
present time. It is obvious that 

UA k A = -N(t)k°. (61) 

A unit vector, tangent to the spatial homogenous slice today t Q , in our position O, can be 
written as 

(k A ) 

e A = (u A ) + / J° = sin 9 cos u> Y[ + sin 9 sin <p Y 2 + cos 9 Y 3 , (62) 
^ >o (u A k A ) o 

where the sky angles 9, ip specify the direction of the incoming photon. Note that we 
introduced on the spatial hypersurface an orthonormal basis {Yi} with g{Yh Yj) = $ij ■ ^ is 
related to the invariant basis {cr a } : 

a 1 = e ~ xl dx 2 , a 2 = e^'dx 3 , a 3 = dx l . (63) 

X 1 = e-* 1 d 2 , X 2 = e~ xl d 3 , X 3 = dx. (64) 
by the following relations 
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Yl 



Y 



-*-X 3 , Y 3 — — — I X 2 Xi 

Jo Jo \ O'o , 



where di = d/dx 1 . The symbols a Q , b Q , f Q denote the scale factors evaluated at present. We 
can parametrise our position without losing the generality demanding (x t ) = .Then we 
get 



Yl 



:<9 2 , Y 2 



1 



<9i, Y 3 = ^-\d 3 - — d 2 

Jo V a ° / 



/Qo VJo 

Finally we can express the present components of the tangent to the null geodesies vector 
from Eg.©, Eq.(fTD 



(k 1 ) = -(k°) J-T sin 6 sin ip , 

V Jo 

(k 2 ) = —(k°) \ — \ sin 9 cosy? — — cos^ ] , 



(65) 



(k 3 ) o = -(k ) o y^cos6. 

Jo 

Our purpose is to find the distortion that the anisotropy of the homogenous spatial 
slices induce, in the propagation of the photons and consequently in z e (8,(p) . It is well 
known that the inner product of a geodesic field k A with a Killing vector field ^ A is constant 
along the geodesies. In our case, there are three Killing vectors tangent to the homogenous 



hypersurfaces |20 



6 = d 2 , 6 = <9 3 , & = d 1 + x 2 d 2 + x 3 d 3 



(66) 



The three constant quantities, along the geodesies, K a = qab £,ak B can be used to determine 
the null vector as 



k 1 = 7~ 1/3 (K 3 - K x x 2 - K 2 x 
k 2 = 7~ 2/3 e 2xl (K ll22 - K 2ll2 ) 
k 3 = 7 ~ 2 / 3 e 2xl (K 2lll - K lll2 ) . 

From Eq.(^) we can compute the three constants 



(67) 
(68) 
(69) 
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A' 



A 1 = ^A 2 + _£«L_ (sin 9 cos ^ - b olo ~ 1 ^ cos 



A 2 = c, 
A3 = 7 2 ^ 3 sin # sin ip . 



Fr^b 



a lo 1 ^ cos 9 + J — — (sin cos <p — & o7o ^ 3 cos 



(70) 
(71) 

(72) 
(73) 



If £ is an affine parameter along the null geodesies then k A 
and o?s 2 = , give using Eq. (|TD|) , and Eqs.( |B7| - |B5D . 



E(9,ip,t a ;t) 



k° 



r 1 " (As-Aix 2 



^E(9, <p,t ;t) 



dx A /dC, . The last relation 



(74) 



A 2 x J 



+ e 2^ 7 -2/3 



(A 2 722 + A 2 7ll - 2AXA27X2 



nl/2 



(75) 



Finally, we find from Eq . ([74]) and Eqs.( |67l - |69l ) that the x l along the geodesies are given by 

- = -7W 1/6 ^y ( A 3 " Aix 2 - A 2 * 3 ) (76) 



dx 
~dt 
dx 2 
~~dT 

dx 3 



-7(*) 



-1/6. 



5(f) 



e 2 * 1 ^ c(f) - A 2 b(t)] 



(77) 
(78) 



This system of equations together with the initial condition x l {t ) = x % = , 7(t D ) = 7o (see 
Eq. |49|) can be integrated numerically. The results of this integration will then be used to 
evaluate the coordinates x l (t e ) of photons, at the emission time t e , and then their redshift. 
From Eqs. (0), (J6T|) , (|74|) we find 



l + z e {9,ip) 



T 



7 ; 1/2 E(9,(p,t ;t e ) 



(79) 



T (e,<p) 

Note that the dependency on the parameters of the finite coordinate transformation is 
manifest in the equation Eq. (|T9"D . In the diagonal form the system of equations describing 
the null geodesies can be obtained from the system Eq.fl3"2p, Eqs.( |76}[77H78] ) effectively, by 
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chopping off the terms containing b, b a . The observed quantity ^f- is related with the above 
expressions through 



(80) 



T (T /T e ) 

where the mean value is understood as an integration over all the sky angles 9, p. The 



quadrupole moment of the temperature field ^ is given by 



a 2 



where 



167r 3 / 2 |(T /T e 



3 (/r+/ 2 2 )+i2 (ii+ii 



T 2 



1/2 



(81) 



h = dip r d9 sin 3 9 sin 2p-^ (9, <p) 
Jo Jo T e 

I 2 = dp r d0 sin 3 9 cos 2p^(9,p) 
Jo Jo T e 

1-2-K PIT J" 

/ dp I dO sin 2 9 cos 9 sin ip-^(6, p) 
Jo Jo T e 



h 
h 



2tt pit f 

dp / d9 sin 2 9 cos 9 cosp— °(9, p) 



o Jo 

2n t-K 



J 5 = dp r d9sin9 (3cos 2 9 -l)^(9,p) 
Jo Jo v ' T e 



and 



(To/T e ) = r dp r d9sm9 ^(0,<p). 
Ait Jo Jo ± e 



(82) 
(83) 
(84) 
(85) 
(86) 



(87) 



The results of the numerical work done, for the general Bianchi V cosmological solution with 
dust, follows: 

• The solution in the off diagonal form has five arbitrary parameters. We determine two 
of the five parameters, choosing values for q Q , Q Q , H Q . Doing this we also determine 
the present value of the shear anisotropy \^j t d = which is small if q D ~ 0, o /2. 
We are left with three free parameters k, d\. These extra three parameters, which 
are not present in the diagonal dust Bianchi solution, studied in 0, determine how 
much the diagonal elements differ inter se as well as with the non diagonal element. 
The choice of different values for these extra parameters select a solution in the space 
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of the g.c.t. equivalent metrics. However these different choices must not change 
the multipole moments since they are true scalar numbers. Indeed the numerical 
calculation exhibits exactly this; the quadrupole moment a^, for the diagonal and 
the non-diagonal form is the same up to next to next order of magnitude. We find 
that the anisotropy and consequently the temperature fluctuations induced from the 
geometry, are small if the shear anisotropy ( U ) and the electric and magnetic scalars 



^ ^uE^^^lj are extremely small. All these three scalars are attaining 
very small values if the value of w is close to zero. 

Compared to the previous results in we found somewhat smaller temperature fluc- 
tuations induced from a dust Bianchi V solution to the CMBR for (^j = 10~ 9 . We 
computed the temperature map and we found that, both for the solution in the diag- 
onal form, analogous to the solution used in 00, and for the solution written in the 
off diagonal form, the temperature allover the sky can be T Q = 2.7277 ±10 -5 when 
(^j = 1CT 9 . As we have already pointed out, this numerical value does not depend 
on the choice of values for the three free parameters k, $\, $\ . 

In solving numerically the system Eq. (|32|) , Eqs.( [?6| - |77| - |78"|) one has to exercise great 
care and demand large accuracy! since otherwise large anisotropics, of non geometric 
origin, can appear as a result of the poor precision. 



The anisotropic distortion (^) rm . = \j ^ it dip £ d6 sin 9 (^) ,of the CMBR is 
smaller for small values of the flatness parameter Q as can be seen below. 



1 However the results in this work, are given for matter content, consisting of dust and radiation. 
2 Large working precision is required, much more than the double precision (16-digits) normally 
used. 
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We found small values of ai coming after integrating all the sky. 
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An interesting outcome of our calculations is that the solution in the off diagonal form 
yields a temperature map all over the sky that is different from the one resulting in the 
diagonal form. The dust diagonal metric, we used, gives a quadrupole like structure 
presented in Fig. 1, while the non diagonal metric gives a quadrupole like structure 
given in Fig 2. Both Figures are drawn for Q a = 0.4 and y^J = 10— . Note that 
the actual observed angles in the sky are shifted from the ones used in calculations 

®obs = 71-0, fobs = TT + (fi ■ 



As far as the focusing effect is concerned, known to appear in open anisotropic models 



21]], the depicted situation inferred from Figure 1 is in complete agreement with the 
findings in |J; the distance between the focusing point and the points of the maximum 
anisotropy is about tt/8. If one wants to reproduce the corresponding Figure in ||, 
from our Figure 1, one must shift properly the axes (the focusing point which is close 
to# = 7r/2, ip = tt/2 should be shifted to 9 = 0, (p = ). Lastly we observe that in 
Figure 2 one gets the same magnitude of the focusing angle. 



IV. CONCLUSIONS 



The exact solution of a Type V Bianchi Universe consisted of an untilted perfect fluid 
with a diffuse matter's equation of state is studied. A simple method for solving the null 
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geodesies is employed, which makes use of some appropriate inner products, constant along 
the geodesies. 

We found that small temperature fluctuations consisted with the observational data can 
be achieved if the parameters of the model are such that give (^j < 1CT 9 . This fact justifies 
our assessment that the model presented, can be a realistic description of our Universe at 
least for the epoch commencing from last scattering and reaching the present time. 

The most important feature of the model studied is that the temperature pattern pro- 
duced by the geometric anisotropies, repeatedly calculated in the literature, estimated after 
solving the null geodesies, is a strongly dependent on the spatial coordinate system used, 
picture. We exhibited how large this dependency is. 
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FIGURES 

FIG. 1. The temperature map for the diagonal solution with Q Q 
temperature values are lighter. 



= 0.4. Regions with higher 



FIG. 2. The temperature map for the non diagonal solution with f2„ = 0.4. 
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